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Abstract 

This paper applies the recently developed theory of discrete nonholonomic mechanics to the 
study of discrete nonholonomic left-invariant dynamics on Lie groups. The theory is illustrated 
with the discrete versions of two classical nonholonomic systems, the Suslov top and the Chap- 
lygin sleigh. The preservation of the reduced energy by the discrete flow is observed and the 
discrete momentum conservation is discussed. 
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1 Introduction 



The theory of variational integrators for Lagrangian and Hamihonian systems originated in j22| . 
125], and It was further developed by a number of authors (see e.g. 0], ^3] ^Kj: 121] > and 
[T7) for a more complete list of references and history). A very important feature of variational 
integrators is the discrete momentum preservation: if the original continuous-time system has a 
symmetry and conserves the momentum map, so does the associated discrete-time mechanical 
system. 

In |13| the theory was extended to the Lagrangian systems with nonholonomic con- 
straints. In particular, it was shown in (2| that the discrete-time nonholonomic system conserves 
the spatial momentum in the case of horizontal symmetry (see |H] for the definition of the hor- 
izontal symmetry). However, the case of horizontal symmetry is not typical in nonholonomic 
mechanics. Apparently, Chaplygin was the first to observe the link between symmetry and 
conservation of the components of momentum relative to the moving frame (see also |28) and ref- 
erences therein). Therefore, it is natural to ask whether the discrete momentum is preserved by 
the discrete-time nonholonomic system associated with a momentum-preserving continuous- time 
system. A closely related property is the existence of an invariant measure. The continuous- time 
nonholonomic systems generically are not measure-preserving (see ^.nd j2Zl for details). The 
next version of this paper will address the measure-preservation property for the discrete-time 
nonholonomic systems. 

The goal of this paper is to study the properties of the numerical variational integrators for a 
nonholonomic mechanical system whose configuration space is a Lie group G. Here we consider 
LL systems, that is, we assume that both the Lagrangian and the constraint distribution are 
invariant with respect to the induced left action of G on TG. 

The paper is organized as follows: Section 2 gives a brief overview of both continuous and 
discrete-time nonholonomic systems. 

In Section 3 we develop the theory of discrete left-invariant nonholonomic systems on Lie 
groups G. The fact that the constraints on G x G are left-invariant enables us to reduce 
the dynamics on a smooth admissible displacement subvariety 5 C G, which is chosen to be 
the exponent of a linear subspace d of the Lie algebra g of G. Under the discrete Legendre 
transformation, S gives rise to a discrete momentum locus U in the coalgebra g*. In contrast 
to continuous nonholonomic systems, the locus is not a linear subspace in g*, but rather a 
nonlinear subvariety. The dynamics is then described by the discrete Euler-Poincare-Suslov 
equations that generate a (generally multivalued) map from U onto itself. 

In Sections 4 and 5 we review the dynamics of the two classical nonholonomic LL systems 
on the Lie groups SO{3) and SE{2), the Suslov problem and the Chaplygin sleigh respectively, 
as well as their multidimensional generalizations. 

In Sections 6 and 7 we construct the discretizations of the above problems as multi-valued 
maps on certain two-dimensional non-orientable subvarieties of S'0(3) and SE{2). It is shown 
that the discrete model retains such a distinct feature of the continuous-time dynamics as the 
existence of heteroclinic trajectories that connect the two one-parameter families of relative 
equilibria of the system. If, for special values of parameters, the continuous-time system is 
momentum/measure preserving, then so is its discrete analog. 

Moreover, it appears that in both discretizations the corresponding reduced constrained 
energy is preserved as well. This conservation law replaces the momentum conservation in 
the general case and seems to be quite unexpected, since generically the discrete variational 
integrators do not preserve the energy and this property does not change in the nonholonomic 
case. 

2 Lagrangian Mechanics with Nonholonomic Constraints 

In this section we briefiy discuss the main concepts of nonholonomic dynamics. For a complete 
exposition see and 0- 
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The Euler Lagrange Equations for Nonholonomic Systems. A nonholonomic 
Lagrangian system is a triple (Q, L, V), where Q is a smooth n-dimensional manifold ealled the 
configuration space, L : TQ — > K is a smooth function called the Lagrangian, and V C TQ is 
a fc-dimensional constraint distribution. Let q = (g^, . . . , q") be local coordinates on Q. In the 
induced coordinates {q,q) on the tangent bundle TQ we write L(q,q). It is assumed that the 
Lagrangian is hyperregular, i.e., the map 

oq 

is invertible (see [TC]'). 

A curve q(t) € Q is said to satisfy the constraints if q{t) £ for all t. The equations 

of motion are given by the following Lagrange-dAlembert principle: The Lagrange-d'Alemhert 
equations of motion for the system are those determined by 

\{q\q')dt^O, (2.1) 

where we choose variations Sq{t) of the curve q{t) that satisfy Sq{a) — Sq{b) = and Sq{t) G 'T>g(^t^ 
for each t where a < t < b. This principle is supplemented by the condition that the curve itself 
satisfies the constraints. Note that we take the variation before imposing the constraints; that 
is, we do not impose the constraints on the family of curves defining the variation. This is well 
known to be important to obtain the correct mechanical equations (see Bloch, Krishnaprasad, 
Marsden, and Murray 3 for a discussion and references). 

Assuming that the constraint distribution is specified by a set oi n ~ k differential forms 

V = {qeTQ\{A^iq),q)=0, j = 1, . . . , s = n - k}, (2.2) 
equation (|2.1|l implies 

d dL dL , ,0/ N „x 

Equations (|2.3|l are called the Euler-Lagrange equations with multipliers. Coupled with l|2.2|l . 
they give a complete description of the dynamics of the system. 

Lemma 2.1. Equations ()2.3f) conserve the energy 

^ = (2.4, 

Proof. Differentiating 12.411 along the flow (|2.3|) . one obtains 

■ / ddL \ /dL .A /dL .\ /dL .. 
^ \dt^qV + W'V W'' 
ddL dL .\ ^ . .\ ^ 

since {A^{q), q) — 0, j = 1, . . . , s. □ 

The Euler^Poincare— Suslov Equations. Let the configuration space be an n-dimensional 
connected Lie group G with local coordinates g. Let g be the Lie algebra of G, that is, the tan- 
gent space TgG at the identity element e G G supplied with an antisymmetric bracket operation 
h-] : X ^ 0. 

Define an LL system on G as a Lagrangian system (G, L, T>) with a left-invariant Lagrangian 
L : TG — !■ M and a left-invariant (generally nonintegrable) distribution T> on the tangent bundle 
TG. 
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The Lagrangian L : TG — )■ R is left-invariant if and only if L{g,g) depends on {g,g) through 
the combination to = g^^g, i-e., there exists a function / : g — > R called the reduced Lagrangian 
such that L{g,g) = l{uj). 

A distribution V C TG is left-invariant if and only if there is a subspace C g such that 
Vg = TLgd C TgG ioY any g g G. Let g* be the dual of the Lie algebra and , j ~ 1, . . . , s, 
be independent elements of g* associated with the subspace 3, i.e., 

f = {ee0l (a^0 = 0, j = l,...,s}- 

Then the left-invariant constraints can be written as 

(a^c^)=0, (2.5) 

where to — g'^^g = TL^-ig is the body velocity operator. 

Define the body momentum p : g* — > M by the formula p = dl/duj. According to |12j . the 
reduced dynamics of an LL system (G, i, T)) is governed by the Euler-Poincare-Suslov equations 

s 

p = ad>-H^A%j (2.6) 

coupled with the constraint equations (|2.5|) . The dynamics of the group variables g is obtained 
by solving the reconstruction equation 

g = TLgLo. (2.7) 

Theorem 2.2. The Euler-Poincare-Suslov equations conserve the reduced constrained en- 
ergy 

Proof. To prove this statement, observe that the reduced energy, {p, lu) — 1{lo), equals the energy 
as the Lagrangian is left-invariant. Since a; G throughout the motion, the reduced constrained 
energy equals the energy along the trajectories of 1)2.6(1 and therefore is preserved. □ 

Let the reduced Lagrangian l{uj) be the quadratic form I = ^{lo,1lo), where I : g ^ g* is a 
symmetric non-singular inertia operator. In this case p = Iw. Then the constraints H2.5|l imply 
that p lies in the subspace 

5^ = {(a^rV) -0, j = l,...,s}cg*. 

It is often convenient to choose a basis ei, . . . ,e„ in the Lie algebra g such that = e"^^"'"^, 
j — 1, . . . ,s. In such a basis, the reduced constrained energy becomes 

i ^ ^u:'u;^ = i ^ r^p.p,. (2.8) 

Here and elsewhere, the quantities J'-* represent the components of the inverse constrained 
inertia operator 

Remark. In the absence of constraints, equations H2.6|l become the Euler-Poincare equations, 
which conserve the spatial momentum J — Ad* p. In the presence of nonholonomic constraints, 
neither the spatial, nor body momentum is conserved generically. However, in some cases the 
body momentum is preserved. (The conditions for the body momentum preservation can be 
seen in EH]). 
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Discrete Mechanical Systems with Nonholonomic Constraints. According to 
a discrete nonholonomic mechanical system on Q is specified by 

(i) a discrete Lagrangian Ld '■ Q x Q R; 

(ii) an (n — s)-dimensional distribution V on TQ; 

(iii) a discrete constraint manifold Vd d Q x Q, which has the same dimension as V and 
satisfies the condition {q, q) E Vd for all q E Q. 

The dynamics is given by the following discrete Lagrange-d'Alembert principle (see j^), 

JV-l , . 

iDiLd{qk,qk+l) + D2Ld{qk-l,qk)j Sqk ^ 0, SqkEVg^, (qt, gt+i) G Pd- 
k=Q ^ ^ 

Here DiLd and -D2-^d denote the partial derivatives of the discrete Lagrangian with respect to 
the first and the second inputs, respectively. 

The discrete constraint manifold is usually specified by the discrete constraint Junctions 

J^,{qk,qk+l)^0, j^l,...,S, (2.9) 

which impose the restriction {qk,qk+i) G on the solution sequence {{qk, qk+i)}- 

Remark. If the discrete Lagrangian is obtained from a continuous one, L{q,q), via a dis- 
cretization mapping : Q x Q ^ TQ defined in a neighborhood of the diagonal of Q x Q, i.e., 
L(j = L o vp, then the variety Vd must be consistent with the continuous distribution V: Vd is 
locally defined by the equations z/^ o v[/ = Q, j = 1, . . . , s. We emphasize that the discretization 
mapping is not unique and hence there are many ways to define the discrete Lagrangian Ld and 
the discrete constraint manifold Vd for a given nonholonomic system {Q,L,V)} 

The dynamics of a discrete nonholonomic system is represented by sequences {{qk.qk^i)} 
that satisfy the discrete Lagrange-d'Alembert equations with multipliers 

s 

DiLd{qk.qk+l) + D2Ld{qk-l,qk) ^^\]A^ {qk). Tj{qk, qk+l) ^ 0- (2.10) 

As in the continuous-time case, these equations are equivalent the discrete Lagrange-d'Alembert 
principle. 

Remark. According to "F, equations H2.10() introduce a well-defined mapping {qk-i,qk) 
(qk, qk+i), if the (n + s) x [n + s) matrix 

/DiD2Ld{qk,qk+l) A^iqk) ••• A' {qk)\ 
D2Ti{qk,qk+l) ••• 

V D2J's{qk,qk+l) ••• y 

is invertible for each {qk,qk+i) in a neighborhood of the diagonal of Q x Q. 

3 Discrete Euler— Poincare— Suslov Equations 

Continuous and Discrete Left-Invariant Lagrangians. Assume that the configu- 
ration space is a Lie group G and denote the local coordinates in G by g. Let the discrete 
Lagrangian Ld G x G G he invariant with respect to the left diagonal action of G on G x G: 

Ld{g9k,ggk+i) = Ld{gk,gk+i) 

^An alternative approach to the discretization of nonholonomic systems based on a modification of canonical 
transformations was proposed in 
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for any g ^ G. 

Define the incremental displacement by the formula Wk — g^^gk+i G G. Since Ld is left- 
invariant, there exists a function Id : (G x G)/G' = G — > R called the reduced discrete Lagrangian 
such that Ld{gk,gk+i) = ld{Wk)- 

According to . for a given continuous left-invariant Lagrangian L{g,g) — l{g~^g) its 
discrete analog Id can be chosen in form 

ld^liilogWk)/h), 

where log : G ^ g is the (local) inverse of the exponential map exp : g G and h £ R+ is the 
given time step. 

For a matrix group G, one can approximate (logWk)/h with 
Wk- I = gk^igk+i - gk), so that 

Ld{gk,gk+i) = I {9k^{9k+i - gk)/h) . (3.1) 

This will be our default choice for the groups SO{n) and SE{n) considered in the next sections. 
Similarly to 0] 115) , we define the discrete body momentum : G x G — > g* by the formula 

Pk = LlD2Ldi9k^i,gk) = L*w^Dld{Wk). 

where I/*^ is the induced left action L*^ : T*G g*. Equivalently, Pk is defined by any of the 
conditions: for any w € g, ^ 



{io.Pk) ^-^Ld{gke''',gk+i) 
ae 



or {uj,Pk)^~^ld{e'^'Wk) 

e=o 



(3.2) 

e=0 



In the unconstrained case, this defines the discrete Legendre transformation C : {gk^Wk) G 
G X G — > {gk,Pk) G G X g*. The mapping £ is uniquely invertible in a neighborhood of the set 
{{g, P) & G X Q* \ P = 0}, but it may fail to be globally invertible. 

In the presence of generic discrete constraints 1)2.911 , the displacement Wk is restricted to the 
admissible displacement subvariety 

Vk = {WkeG\ T,{gk,gkWk) = 0, j = 1, . . . , s} 

As a result, the discrete momentum Pk is restricted to an (n— s)-dimensional subvariety C g*, 
the image of C{gk, Vk) in g*. In case of generic discrete constraints, for different k the subvarieties 
Uk are different. 

Discrete Left-Invariant Constraints. If the continuous constraint distribution V is 
left-invariant, it is natural to require that the discrete constraint manifold T^d is invariant with 
respect to the left diagonal action of G on G x G, that is, 

^]{9 9k,99k+i) = ^j{qk,qk+i) for any g E G, j = l,...,s. 

This implies that there exist functions /j : G ^ M, j = 1, . . . , s, such that 

Tj{qk,qk+l) = fj{Wk). 

Consequently, Vd C G x G is completely defined by the admissible displacement subvariety 

S = {h{W) = 0, MW) = 0}cG, 



^The definition of the discrete momentum 13.211 accepted in this paper computes Pk as a function of Wk whereas 
the standard definition used in many publications, 

{u;,Pk} = -^Ld {gk-\,gke'^' 

makes pk a function of Wk-i- 
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namely Vd = {g, gh}, g eG,h eS. 

The submanifold S should pass through the identity element / in G, and the tangent space 
TSi at the identity should be "horizontal", i.e., it should coincide with the linear subspace 
c) C generating the left-invariant distribution on TG. 

The second property suggests that S — {W G G \ logM^ € O}. 

Equivalently, S can be chosen a union of all one-parameter subgroups G^^ generated by all 
admissible vectors ci) G 0. In other words, one can set S — expO. However, in case of generic 0, 
the set expO is not a subvariety of G. 

In this paper we concentrate on the important case when G contains a subgroup H generated 
by subalgebra () C g such that there is a decomposition g = i) ® and (g, f)) forms a symmetric 
pair, that is 

[f),f)]cf), [0,0] Ct], [f),f>]cf.. (3.3) 

Proposition 3.1. (^01) Under the condition the set S — exp is a smooth submanifold 

of G, which is either homeomorphic to the symmetric space G/H or is a factor of G/H by a 
finite group action. 

Under conditions H3.3|) the set expO is known as the Cartan model of the symmetric space 
G/H. 

Notice that the tangent bundle T 5 is not a subset of the left-invariant distribution D C TG, 
since the latter is not integrable. 

Under the Legendre transformation C, the discrete momentum P^. is confined to the subva- 
riety 

U = {peg* \p = LMW),W eSjdg*, 

which now docs not depend on fc. It appears that in the examples considered below the map 
S t-^ U is uniquely invertible almost everywhere on U. 

Discrete Euler— Poincere— Suslov Equations. Assume that the discrete Lagrangian 
: G X G ^ R, the discrete constraint distribution Vd, and the constraint distribution D are 

left-invariant with respect to the left action of G on G x G and TG, respectively. 
Define the action .sum and the reduced action sum by the formulae 

N-l N-1 

Sd^^ Ld{gk,gk+i) and = ^ ld{Wk), 

k=0 k=0 

respectively and rewrite the nonholonomic constraints H2.2|l as a set of vanishing one-forms 

Following [H| , consider variation of Sd assuming that the variations Sg^ satisfy the conditions 
{gk)Sgk = 0, j = 1, . . . ,s and Sgo = SgN = 0. 

For the left-invariant constraints given by H2.5|l the admissible discrete variations are those 
Sgk G TGg^ that satisfy the conditions 

{a^,g^'Sgk) =0, j^l,...,s, k^l,...,N-l. (3.4) 

The following theorem extends the result of to the nonholonomic setting. 

Theorem 3.2. Let Ld : G x G ^ M. be a left-invariant Lagrangian, Id : G ^ M be the 
reduced Lagrangian, andV and Vd be the compatible constraint distributions on TQ and Q x Q, 
respectively. Then, following statements are equivalent: 

(i) The sequence {{gk, gk+i)}^^i^ is a critical point of the action sum Sd ■ G^^^ —* R for 
arbitrary constrained variations. 

(ii) The sequence {{gk, gk+i)}^^Q satisfies the discrete Euler-Lagrange equations with multi- 
pliers Y2.10\) with q replaced by g, that is, 

s 

DiLd{gk,gk+i) + D2Ld{gk-i,gk) ^^K^jidk) (3.5) 

i=i 
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which are coupled with the discrete constraint equations J^j{gk, gk+i) = 0. 
(Hi) The sequence {Wk}^SQ is a critical point of the reduced action sum Sd ■ G^'^^ R with 
respect to variations SWk, induced by the constrained variations Sgk, and given by 



9kli^9k+i - Ad^-i ^Sgk ■ (3.6) 



(iv) The sequence {Wk}i^^Q satisfies the equations 

s 

QWk-i)TLw,_, ~ l'^{Wk)TLwMw-^ = E (3.7) 

coupled with the discrete constraint equations 

/,(W^fe) = 0, j = fc = l,...,iV-l. 

Proof of the theorem is given in the end of the section. 

We now rewrite (|3.7|) in the form of discrete momentum equations. For any w e g, 

UWu-i)TLw,_,u: = {UWk-i),TLw,_,u:) = {TL*^^J^{Wk-i)M 

and similarly 

l'^{Wu)TLwMw-- ^ = {l'AWk),TLwMw, ^) = (Ad;^-i T L\yJ'^{Wu) , u:) . 

Therefore, in view of the definition of the discrete momentum (13. 2|) . ^A.7\ becomes discrete 
Euler-Poincare Suslov equations 

s 

where Wk is restricted to S and pk <^U <Z Q* . 

The above equations extend the discrete Euler-Poincare equations obtained in |^ ^] to the 
case when the discrete left-invariant constraints are present. Thus, they represent a discrete 
analog of H2.6|l and define a map B : U ^ U : Pk ~* Pk+i, which is generally multi- valued. 
Given Pk, one evaluates Pk+i by 

1. Finding Wk by inverting the Legendre transformation; 

2. Calculating Pk = Adjy^ Pk; 

3. Choosing Pk+i as one of the points of intersection of the {n — s)-dimensional subvariety U 

with the linear space span(ai, . . . , as) passing through Pk- 

If the map is multivalued, one needs to make a choice of a branch of B. One natural way of 
doing this is to start from a value of Pk whose norm is small and to select Pk+i of the smallest 
norm. 

Proof of Theorem \3./l\ We first prove the equivalence of (i) and (ii) following ;6 . Recall that the 
variations 5gk vanish at /c = and k = N . Computing the first variation of the discrete action 
sum Sd, we obtain 

N-l 

SSd = E ^d.{9k,9k+i) 

k=0 

N-1 N-1 



^ DiLd{gk,9k+i)S9k + E D2Ld{gk,9k+i)5gk 

k=0 k=Q 
N-1 N-1 

^ DiLd{gk,9k+i)S9k + E D2Ld{gk^i,gk)Sgk 

k=l k=l 

N-1 

E {DiLdi9k, gk+i) + D2Ld{9k-i, gk)) Sgk ■ 



k=l 
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Here the variations Sgk are not independent and satisfy the conditions Aj{gk)5gk = 0. Therefore, 
SSd = if and only if 1231) is fulfilled. 

Next, we prove that (i) is equivalent to (iii). Notice that Ld — Id ° t^, where tt : G x G ^ 
(G X G)/G = G is given by {gk,gk+i) 9l^9k+i- Therefore 



5sd = 5Sd 



The variation 5Wk is computed to be 



5Wk = S{g^^gk+i) = g^^Sgk+i + Sg^^gt+i 
= 9k^Sgk+i - 9k^5gkgk^gk+i 

= i9k^9k+i){gkliSgk+i) - {9k^9k+i){9k^gk+iy^{gk^6gk){gk^gk+i) 



9k 9k+i 



9kli59k+i - Ad(^_,^^^^)-i (5, 



which yields 

To prove the equivalence of (iii) and (iv), we use (|3.t)|l to compute 



N-1 



N-1 



6sd = Sj2 IdiWk) = J2 ^^diWk) = ld{Wk)SWk 

k=0 fc=0 k=0 

= J2 ^'di9k^9k+i)9k^9k+i [g^l^Sgk+i - Ad, {g^^Sgk) 

k=0 
N 

^Y^'d i9k-i9k) {9k-i9k) {9k^5gk) 

Y (.9k'9k+i) ((.9fe \9fc+i) Ad^^_i^^^^^-i {g^^Sg, 



k=l 



k=0 



N-1 



E [^'d i9k\9k) TL^-^^^^ 

I'd {g^'gk+i)TL^-.^^^M^^^,^^^^y 



fe=i 



{9k <^9k)- 

Since the variations Sgk satisfy the conditions ()3.4|) . Ssd = if and only if item (iv) holds. 



4 The Suslov Problem and its Multidimensional General- 
izations 

The most natural example of LL systems is the nonholonomic Suslov problem, which describes 
the motion of a rigid body about a fixed point under the action of the following nonholonomic 
constraint: the projection of the angular velocity vector lo € to a certain fixed in the body 
unit vector 7 equals zero: 

(^,t) = 0. (4.1) 

The configuration space of the problem is the group 50(3). Under the identification of Lie 
algebras (M"^, x) and (so(3), [•,•]), oj and 7 correspond to elements of so(3) and the coalgebra 
so* (3) respectively. 

Let I : M'^ i-^ K."^ be the inertia tensor of the body. Then the Lagrangian equals L = 
^(t?,Ia;) and the momentum p is represented by the vector M = {Mi, M2, M^)'^ = Tuj. The 
left action of the group 50(3) on r50(3) leaves the kinetic energy of the body and the constraint 
(|4.1|l invariant. 
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For the Suslov problem the Euler-Poincare-Suslov equations (|2.()|) on so(3) become 



^(ITZf) = IIj X + At*, (4.2) 

where x denotes the vector product in R'^ and A is the Lagrange muhipher. Differentiating 
<|4.1(l . we find 

A = -{ITS X T?,ri7)/(7',ri^). 
Therefore, H4.2|l can be represented as 

|(i^) = ^:^^rWx((i7?x^)xt), 

which, in view of H4.1|l . is equivalent to 

— (W) = (IT? , ^)Tj X I-i^. (4.3) 
dt 

The Suslov system possesses the energy integral 

(7j , ITj ) = (M, r^M) =h, const (4.4) 

and, as seen from H4.3|l . it has a line of equilibria positions 

i? = {(i:?,^) = o}n{(i7i;,^) = o}. 

Note that in the principal basis, where I — diag (Ii, I2, 13), the system has the integral given by 
degenerate quadratic form 

(I27I + II372 ^137172 ^I27l73 \ 
-I37172 II17I + 137? -I17273 , (4.5) 
-I27173 -II17273 I172+I271/ 

which coincides with the restriction of (|4.4(l onto the constraint plane (M,I^^^) = 0. 

In the basis where only one of the components of ^ is nonzero, say ^ — (0, 0, 1)^, and 
the inertia tensor is unbalanced, the integral H4.4|l can be replaced by the reduced constrained 
energy integral 

h2Mf - 2I12M1M2 + Mil (4.6) 

The dynamics of the two independent momentum components. Mi and M2, is illustrated in 
the Figure Erl Because of the conservation law H4.6|l . the trajectories are the elliptic arches that 
form the heteroclinic connections between the asymptotically stable (filled dots) and unstable 
(empty dots) equilibria. 

As a result, the motion of the rigid body is the asymptotic evolution from a permanent 
rotation about an axis fixed in the body frame to a permanent rotation about the same axis 
and with the same angular velocity, but in the opposite direction. Note that in space the axes 
of the limit permanent rotations are different. 

The Suslov problem admits some natural multidimensional generalizations studied in |71 
13 I26| . The configuration space of an n-dimensional rigid body with a fixed point is the Lie 
group SO{n). For a path R{t) e SO{n), the angular velocity of the body is defined as the 
left-trivialization uj{t) — ■ g{t) e so{n). 

The left-invariant metric on SO{n) is given by non-degenerate inertia operator I : so{n) 
so{n). Then the Lagrangian of the free motion of the body reads 

L^^{lcj,cj), (4.7) 
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Figure 4.1: The Momentum Dynamics for the Suslov Problem. 



where now (•, •) denotes the KiUing metric on so{n), {X, Y) = — ^tr {XY), X,Y E so{n). For a 
"physical" rigid body, luj has the form Joj + ojJ, where J is a symmetric n x n matrix called 
mass tensor (see 0). 

Let ei , . . . , e„ be the orthogonal frame of unit vectors fixed in the body. What form may have 
a multi-dimensional analog of the condition l|4.1(l ? To answer this question, note that, instead of 
rotations about an axis in the classical mechanics, in the n-dimensional case we have infinitesimal 
rotations in the two-dimensional planes spanned by the basis vectors Ci, ej, i, j — 1, . . . ,n. 

Suppose, without loss of generality, that 7* — (0,0,1) in H4.1|) . Then this condition can 
be redefined as follows: only infinitesimal rotations in planes (61,63) and (62,63) are allowed. 
Hence, it is natural to define the n-dimensional analog of Suslov's condition in the following way: 
only infinitesimal rotations in the planes (61, 6„), . . . , (6„_i,6„) (i.e., in the planes containing 
the vector 6„) are allowed. Thus, in the above basis, the angular velocity matrix in the body 
must have the form 

/ ... com \ 







This implies the constraints 

(w, 6i A Cj) = [ci, uJCj) = 0, 1 < i < i < n — 1. 



(4.8) 



(4.9) 



As a result, the multidimensional Suslov problem is described by the EPS equations on the Lie 
algebra so{n) 

Ap,6pA6„ (4.10) 

l<p<Q<n— 1 



dt 



where the multipliers Xpq can be found by differentiating the constraints (|4.9|l . 

Integrability of the system (|4.10|) . H4.9|) was proved, and its geometric properties were studied 
in [7], whereas the reconstructed motion on the group SO{n) was described in |26j . 



5 Chaplygin Sleigh 

Another example of a mechanical system governed by the Euler-Poincare-Suslov equations is 
the so-called Chaplygin sleigh, the system introduced and studied in 1911 by Chaplygin jSj (the 
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work had been actually finished in 1906, see also 

The sleigh is a rigid body moving on a horizontal plane supported at three points, two 
of which slide freely without friction while the third is a knife edge which allows no motion 
orthogonal to its direction, as shown in Figure 1531 




Figure 5.1: The Chaplygin Sleigh 



The configuration space of this dynamical system is the group of Euclidean motions of the 
two-dimensional plane M^, SE(2), which we parameterize with coordinates {9,x,y). As the 
figure indicates, 9 and (x, y) are the angular orientation of the blade and position of the contact 
point of the blade on the plane, respectively. 

The Lagrangian and Constraint in the Body Frame. Introduce a coordinate sys- 
tem called the body frame by placing the origin at the contact point and choosing the first 
coordinate axis in the direction of the knife edge. Denote the angular velocity of the body by 
u> = 9, and the components of the linear velocity of the contact point relative to the body frame 
by vi,V2- The set (w, ^^i, ■^2) is regarded as an element of the Lie algebra se(2). 

The position of the center of mass is specified by the coordinates (a, h) relative to the body 
frame (we not assume here that the center of mass lie along the blade direction as in some 
models). We will see that a is crucial to qualitative behavior of the system while h is irrelevant. 

The Lagrangian equals the kinetic energy of the body, which is a sum of the kinetic energy of 
the center of mass and the kinetic energy due to the rotation of the body. Let m and J denote 
the mass and moment of inertia of the sleigh relative to the contact point. The position of the 
center of mass relative to the fixed (inertial) frame is 

{x + a cos 9 — 6sin0, y + asinfl + bcos9). 

Thus, the kinetic energy of the center of mass has the form 

■2 I ■2 I /„2 I i,2\q2 



X 



+ y + (a^ + r)6'^ + 26* (a(-.T sin 6* + y cos 6*) - b{x cos9 + y sin9)) 



or, using the body components of the angular and linear velocity, 

m [(a^ -I- +vl+vl- 2buvl + 2aijjV2] ■ 

As a result, the (reduced) Lagrangian is 

I = i [(J -f m(a^ + b'^))uj'^ + m{vf + v^) - 2mbujvl + 2maujV2\ . (5.1) 
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Next, the constraint written relative to the body frame is = 0. Both the Lagrangian and 
constraint arc invariant with respect to the left action of SE{2) on TSE{2) as they depend on 
(g, g) through the combination 

n = g-^g. (5.2) 

The Dynamics of Chaplygin Sleigh. In view of H5.1|) . the components of the body 
momentum are 

dl 

Pe = = (J + m(a^ + h'^))uj + 2m(av2 — bvi), 
ouj 

dl , , . dl , , 

Pi = ^T— = m(vi - buj), P2 = t;— = m[v2 + auj). 
ovi 0V2 

The reduced dynamics of the Chaplygin sleigh is governed by the equations 

Pe=PlV2-P2Vl, Pl=P2UJ, p2^-piuj + X, (5.3) 

which are the Euler-Poincare-Suslov equations H2.6|l on the algebra se(2) coupled with the 
constraint V2 — 0. Eliminating the variables fl and the Lagrange multiplier A from (|5.3|l . one 
obtains the reduced dynamics of the Chaplygin sleigh in the form of the momentum equation 



(J + ma^)^ 



(pe + bpi) [mbpe + ( J + m(a^ + 6^)) pi) 



ma / , , n2 
(J + ma'^Y 



(5.4) 



which has the constrained energy integral 

mpl + 2bmp0pi + {J + m{a^ + b'^))pl. (5.5) 

In the case 6 = equations 1)5. 4|l become 

apePi . mapl 
Pe = — 5- J, Pi = 7^- jT^- 5.6) 

We emphasize that the phase portrait of l|5.4|l is identical to that in the Suslov problem. 
Indeed, if a = 0, the nonholonomic momentum {p0,pi) is conserved. Therefore, the body angular 
velocity w and the component of the body linear velocity along the blade vi are constants. The 
evolution of the configuration variables {0, x, y) is determined from the reconstruction equation 
H5.2|l . which reads 

6 — uj, i cos 6* + y sin ^? = ui , — a;sin6' + ycosfl = 0. (5.7) 
The solutions of (|5.7() are 

9 — Oo + ujt, X — xq -\ — - sin(6'o + cot), y — yo cos(0o + ^0 if cj ^ 

and 

6 ~ 9q, X = Xq + vi COS 00 1, y = yo + sin 9q if u) = 0. 

Therefore, the contact point of the blade and the plane generically moves along a circle at a 
uniform rate. 

If a ^ 0, the dynamics H5.4|l is integrable as the reduced energy is conserved. The trajectories 
of 1)5. 4|l are either equilibria situated on the line pg + bpi = 0, or elliptic arches.^ The equilibria 
located in the upper half plane are asymptotically stable (filled dots in Figure whereas the 
equilibria in the lower half plane are unstable (empty dots in Figure [^J. The elliptic arches 
form heteroclinic connections between the pairs of equilibria as shown in Figure 15.21 
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(a) a / 0, 6 / 



(b) a / 0, & = 



Figure 5.2: Momentum dynamics. 




(a)a/0,b/0 (b) a / 0, 6 = 



Figure 5.3: Generic trajectory of the blade. 

A generic trajectory of the contact point of the blade and the plane has a cusp point (see 
Figure At the cusp, the speed of the contact point, \vi\, momentarily vanishes as the 

momentum trajectory intersects the line mbpi + ( J + m(a^ + 6^))p2 = 0. 

Since the group SE{2) is a "non-compact" version of the group S0{3), the dynamics of the 
Chaplygin sleigh can be interpreted as a "non-compact limit" of the dynamics of the Suslov 
problem. Recall that the any non-equilibrium trajectory of the Suslov top has a steady-state ro- 
tation as its asymptotic dynamics. In a similar manner, a non-equilibrium state of the Chaplygin 
sleigh asymptotically approaches a uniform straight- line motions as t —^ ±oo. 

The shape of the generic trajectory of the contact point is predetermined by the inertia of 
the body and the position of the center of mass relative to the blade, and is independent of 
the initial conditions. While the dynamics of the group variables (9, x, y) cannot be explicitly 
written, it is possible to compute the angle between the asymptotic directions of the dynamics 
of the contact point. See [S] and for details. 

Multidimensional Chaplygin Sleigh. We now briefly discuss the generahzed Chap- 
lygin sleigh, which is an n-dimensional rigid body moving in in the presence of certain 

^Tiiis follows from matching the trajectories and the level curves of the reduced energy, which is a positive-definite 
quadratic form. 
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nonholonomic constraints. 

The configuration space of this dynamical system is the group SE(n), which has the structure 
of a semidirect product, SE{n) — 50(n)(S)R", so the group elements are written as {R,x), 
where R G SO{n) is the orthogonal rotation matrix of the body and x S R" is the position 
vector of its origin A. It is often convenient to represent the elements of SE{n) by means of 
(n + 1) X (n + 1) matrices of the form 

9{R.x) ={^} \ 

and the group operations for SE(n) correspond to operations with the matrices: the product 
of two such matrices corresponds to the superposition of two Euclidean motions represented by 
these matrices and the inverse matrix correspond to the inverse Euclidean motion. 

The Lie algebra se{n) of the group SE{n) is the semidirect product so(n)(S)M" and it is 
isomorphic to the set of {n + 1) x (n + 1) matrices 

The elements of se{n) are written as The Lie bracket [771,772] in se(n) is 771772 — ^72771, 

which yields 

(6,^2)] = ([6,6],ClW2 ~ ^2Vl). 

For a trajectory g{t) C SE{n), the body velocity operator is defined as the left-trivialization 
^(t) — g^^g{t) E se{n). In this case u! ~ R^^R{t) and v — R^^i{t) are respectively the angular 
velocity matrix and the vector of linear velocity of A in the body frame. 

As in the classical case, we suppose that the center of mass C of the body does not coincide 
with the origin A of the body frame. Let (oi, . . . , a„)"^ be constant position vector of C in this 
frame and, as above, J — diag(Ji, . . . , J„) be its mass tensor. Then the Lagrangian is 

L = -ifr (CJS^) = -^tr {lj{J + ma<» a)uj) +m{v,Lua) + (5.8) 



/I ai\ 

I = S diag (Ji, . . . , Jn, m) S'^ , S — 

1 a„ 

\0 ... ij 



e SE{n), 



where S describes the position of the center of mass C relative to the body frame. 

The body momentum is an element of the dual space se*{n) and it is given by the pair 

P = {M,p) e se*{n), M G so(n), p G M", 
dL dL 

M.,j ^ , p., = — , 7, j = 1, . . . ,7^. 

aujij ovi 

Straightforward evaluation leads to the formulae 

A/ = ( J + 777a ® a)uu + uu{J + ma ® a) + m{v ® a — a®v) G so{n), 

(5.9) 

p^m{v + Loa) 

Here M is the angular momentum of the body with respect to its center of mass C and p is the 
linear momentum of C as a point with mass 777. 



Left-invariant constraints on SE(n). There are numerous ways to introduce nonholo- 
nomic constraints for the generalized Chaplygin sleigh. For example, one can require that the 
velocity of the reference point is restricted to a A:-dimensional linear subspace fixed in the body. 
For 71 = 3, such constraints were studied in [l^j and (27).) 
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Another natural choice is to define the constraint subspace c) € se{n) to be the set of matrices 
of the form 



/ Wi2 

~UJi2 



V 

In the particular case n = 2 we have 



-UJin 





Win Vi\ 






0/ 



(5.10) 



—u) vi^ 
S ^ \ uj 
,0 



(5.11) 



6 Discrete Suslov System on SO{n) 

Now we apply the the discrete Euler-Poincare- Suslov equations H3.8|l to construct a discretiza- 
tion of the Suslov problem. Let Rk G SO{n) be the orthogonal rotation matrix describing the 
k-th position of n-dimensional top. 

Introduce the finite rotation matrix ilk = Rk^k+i, analog of the angular velocity oj in the 
body. Note that in the continuous limit, when Rk+i = Rk + sR, e << 1, one has 

Qk =^1 + R'^R^I + euj, (6.1) 

Define the left-invariant discrete Lagrangian on SO{n) x SO{n) by substituting oj in H4.7fl 
by R^{Rk+i — Rk) = — I- Using the property R^Rk — I, we get 

^d(f^fc) = ^tr(rjfcJ), and ^^(i?^, i?fc+i) = ^ tr(i?fc Ji?^^ J. 

Then, following the definition (|3.2() . the body angular momentum Mk G so*{n) has the form 

Mk — Rk ^Rk+iJ — JRk+iRk = ^kJ — J^k^ (6.2) 

which in the above limit transforms to Juj + to J, the standard relation between the angular 
velocity and momentum. The expressions for Ld, Mk were originally introduced in |18| . 

Remark. In the classical case n = 3 one can parameterize Rk in terms of the Euler angles 
6k,ilJk,4'k, as coordinates on 5*0(3) (see e.g., pS]!. 

cos (pk cos ipk — COS Ok sin (pk sin ^pk — cos 4>k sin ipk — cos 0k sin (j)k cos ipk 
sin (j>k cos Tpk + cos 9k cos (fik sin ipk — sin (j)k sin ipk + cos 9k cos (fik cos ipk 
sin 9k sin tpk sin 9k cos ipk 



sin 6'fc sin 
— sin 9k cos 
cos 9 k 
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Substituting these ones and analogous expressions for Rk+i into the discrete Lagrangian Ld{Rk, Rk+i), 
we obtain 



cos6k cos 6*^+1 [1 + cos(A(/)fe) cos{ipk + tpk+i)] 



+ cos{Tl>k +'0/c+i) sm9k sin6'fe+i + cos(A0fe)[sin6'fe sin6'fe 



+1 



cos(-0fe + -0/c+i)] + -[cos6'fe+i - cos 6lfc] sin(A0fc) sin(i/n; + Vfe+i) 



cos-dk cos?9/i;+i +cos(A(/)fe) cos(A?/'fe) — cosi^fe cosz9fe+i cos(A(/ifc) cos(A'(/'fe) 
+ cos(A0fe) sini?fc sini?fc+i — cos(A-0fc) sini^^ siniJfc+i 
cosi^fc sin((?!)fc + (/)fc+i) sin(A?/'fe) - cosi?fc+i sin(A(?!)fc) sin(V'fc + ■>pk+i) 



A. 



cos(A(/)fc) cos(A'(/;fc) +cos6k cos9k+i cos{A<j)k) cos(A-0fe) 



1 

2 

cos6'fe cos6'fe+i + sin6lfe sin6lfc+i(cos(AV'fe) - cos(A(^fc)) 



(cos6'fc + cos6'fe+i)sin(A(/)fc) sin(AV'fc) 



A, 



(6.3) 



where A^fe = 61^+1 - Ok, A0fc = 4>k+i - AT/jfe = V'fe+i - V'fe, and 
^1 = J2 + J3, ^2 = Ji + J3, ^3 Ji + 

are the principal moments of inertia of the rigid body. 
In the continuous limit, setting in 



'(5t, (l)k+i - (t)k ^ (pSt, il^k+i - ipk = St, St«l 



(6.4) 



then expanding in 6t and dividing by (St)'^, up to an additive constant and terms of order St, 
one obtains the well-known expression for the kinetic energy of the top (see, e.g., |25j) 



T = i(0sin6'sin'i/; + ^cosV')^Ai 

+ i(0sin6'cos^ - 0sini/')^y42 + ^(V^ + cos 6*)^ A3, 



(6.5) 



where the expressions in brackets represent components of the angular velocity vector in the 
frame attached to the body. 

Notice that the discrete Lagrangian (|6.3|l does not coincide with the "straightforward" dis- 
cretization of H6.5|l obtained with a direct replacement of the angular velocities by the angular 
differences according to (|6.4ll . 



Discrete constraints on SO{n). Following the approach described in Section 2, we im- 
pose discrete left-invariant constraints on SO{n) x SO{n) in the form of restrictions on finite 
rotations 0/j G SO{n). In accordance with the continuous constraints H4.9|) . we assume that 
admissible rotations must be exponents of the vectors of the linear space 

d = span{ei A e„, . . . , e„-i A e„} C so{n). 

Lemma 6.1. 1). In the basis ei, . . . ,e„, the admissible rotation matrices have the structure 

{^k)ij = {^k)ji, {^k)in = -ii^k)ni, i<ij <n-l, (6.6) 

that is, they are anti- symmetric in its last row and column and symmetric in the other 
part. 
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2). The admissible displacement subvariety S — expO is homeomorphic to the projective space 
P"^^ = S"~^/'Z?. In the same basis, the components of ilk «re parameterized by points of 
the unit sphere S"^^ — {zq + zf + ■ ■ ■ + z^_i = 1} in the form 

{^k)ij — 5ij — 2ZiZj, {^k)in = -(f^fe)m = 2zoZi, {flk)nn = 2Zq — 1, (6.7) 

l<i,j<n— 1. 

Note that in the continuous limit described by (|6.1() . conditions (|6.6|l yield the constraints 
(|4.8|l on so{n). 

Proof of Lemma 16.11 1). Any vector of C so{n) can be represented in the form 0uAe„, where 9 
is a nonzero constant and u = (ui, . . . , u„-i, 0)"^ is a unit vector in M"^^ = span (ei, . . . , e„_i). 
The odd powers of On A e„ are skew-symmetric and have zero left- upper (n — 1) x (n — 1) part, 
whereas the even powers are symmetric and have zero last row and last column. Hence, the 
exponent of 9u A e„ must be of the form (|6.6|) . 

2). The operator TZe^u — exp(6'u A e„) C SO{n) describes rotation in the 2-plane spanned 
by u, e„ by the angle 6. Then we get 

TZe^u u ~ cosd ■ u — sin6' • e„, 

T^s,u Sj — Cj — (cj, u)u -I- (ej, u)(cos6' • u — sinf? • e„), 1 < j < " — 1, 
T^B.u Cn = cos 9 ■ Cn + sin ■ u. 

The latter n vectors form columns of the matrix TZe^u- Setting in the above formulas 

Zi = sin6'/2 Ui, l<«,i<fT- — 1, 2o = cos0/2 (6-8) 

and identifying TZg,u with flk we arrive at expressions (|6.7|) . 

Since sin(6'/2) = - sin(27r - e)/2 and cos(6'/2) = - cos(27r - 0)/2, from we conclude 

that opposite points on 5"^^ correspond to the same admissible rotation 71$ ^u- Finally, there 
is a bijection between S — expO and P"~i = S^~^ /l? . The lemma is proved. 

Note that (|6.6(l imply left-invariant constraints on SO{n) x SO{n) in the form 
tr(i?fc Cj A e„ Rk+i - R^+i Cj A e„ Rk) = , j = 1, . . . , n - 1. 



Rotations about an axis. In the classical case n = 3 the conditions (|6.6|l say that flk is 
a finite rotation about an axis lying in the plane (ei, 62), while expressions (|6.7|l imply that the 
rotation axis is directed along vector p — (22, — zi, 0)"^ G R^. 

Indeed, the group 5*0(3) can be regarded as covered twice by the unit sphere = {q^ + 
91+92+93 ~ l}j where go, ■ • ■ , 93 are the Euler-Rodriguez parameters such that any rotation 
matrix W G 5*0(3) can be represented in form (see, e.g., ^E^) 

fqo +11- 12-13 2((7i(72 + 9390) -2(9193 - 929o) \ 
W = 2(gig2 - 9390) 1o+4-lf~ ll -2(9293 + 9o9i) • (6.9) 
\ -2(9193 + 929o) -2(9293-9091) 1I+ ll~ I1- I2) 

The operator W describes a finite rotation in M?' about the vector e = (91, 92, 93)^ by the angle 
9 such that 90 = cos 0/2. 

Setting in (|6.9I) {^k)i2 = [^k)2i imphes 93 = 0, hence is a rotation about an axis lying 
in the plane (ei, 62). In this case admissible operators $7 G 5 C 5*0(3) have the form 

(2{ql+ql)-l 2gig2 2go92 \ 

2gig2 2{ql+ql)-l -2go9i , (6.10) 

\ -29 092 2go9i '^ll-'^J 

^Here and below, to simplify notation, we omit the discrete time index k at the components of z. 
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which, under the substitution qi = —Z2, q2 = zi, qq = zq, coincides with the parameterization 
(|6.7() . As a result, the variety of such matrices is the real projective plane RP^ = 

We emphasize that, in general, the k-th position of the body Rk = ^k-i • • • f^o is not a 
rotation in the plane (ei, 62). 

Discrete momentum locus 14 C so* (3). In contrast to the continuous case, the discrete 
momentum Mk does not lie in a linear subspace in the coalgebra so* (3), but on a nonlinear 
algebraic variety U C so* (3) defined by the relation (|6.2() and the conditions H6.1()|l . 

If in the frame 61,62,63 the tensor J is diagonal, J = diag{Ji, J2, Js)^ , then the angular 
momentum vector M — (Mi = — M23, M2 = M13, M3 — —M12Y' has the form 

M = 2((J2 + J3)go9i, (^1 + ^3)9092, (Ji - J2)qiq2 f 

(as above, to avoid tedious notation we omit the discrete time index at the components of q). 
Here and below, without loss of generality, we always assume qo > 0. As a result, U coincides 
with the Steiner Roman surface in R"^ given by the quartic equation 



{J2 + JsK-h + J3) ' " (J2 + J3)(Ji- J2) 

- Af| Mi - 2 Afi M2 Ma = (6.11) 



J2 + J3 



(Ji + JsK-h - J2) 
(see, e.g., IS1I201)- 

In general case, when J is not diagonal in this frame, one has the parameterization 

_^ / {J22 + J3z)qoqi - Ji2qQq2 -{Jiaqi + J23q2)q2\ 
M = 2 \ (Jii + J33)qoq2 - Juqoqi +{Ji3qi + J23q2)qi ■ (6-12) 
\(Jii - J22)gi92 - Ji2(<?i - g|) -(Jiagi + J23'?2)<?o/ 

One can show that the components of M satisfy an algebraic equation of degree 4, which 
generalizes (|6.11|l and which we do not write here. The corresponding algebraic surface lA 
in M?" — {Ml, M2, Ms) has pinch points and self- intersections. One can also show that if the 
quadratic form ( J22 + ^33)'Zi — 2Ji29i(j'2 + ('/33 + '/ii)92 is positive-definite, then any pair (Mi, M2) 
has at most two real inverse images on U. 

An example of such a surface for an unbalanced inertia tensor and its circular section for 
0.4 < qo < 0.8 are given in Figures I^TI 16 . 21 respectivelv. 



Discrete EPS equations on so*(3). In the considered case G — 50(3), the discrete 
momentum equation with multipliers l|3.8|l takes the form 

/O 1 0\ 

Mk+i = nlMknk + Xk [-1 00, Mk^rikJ-jnl, (6.13) 
\ 00/ 

where the components of ilk are subject to constraints (|6.6|l . 

This provides a discrete analog of the Suslov system H4.10|l on so* (3) and defines a map 
U ^ U or, in view of expressions H6.10|l . (|6.12|) . a map 

B : MP^ RP'^ : {qi,q2,qo) {qi,q2,qo), which is generally multi- valued. 
To describe the latter map in details, we note that in (|5.13|) 

fi^Mfefife = jnk - nlj = nlMl 

((J22 + ^33)9091 - Jl2qoq2 +(Jl39l + J2392)g2\ 

{Jii + J33)qoq2 - Juqoqi -(J1391 + -^2392)91 , (6.14) 
-(Ju - ^22)91 92 + Ji2{ql - qj) -{Juqi + -^2392)90/ 
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Figure 6.1: The Momentum Surface U. 



where 'uj denotes vector representation of element uo of so(3). Comparing this with (|6.12(l . we 
find that H6.13|l can be written in form 



k+l 



Mk+4: \ -{Jnqi + ^23'72)'7i 

-(Jii - J22)qiq2 + Ji2{qf - ql) + A*.. 



(6.15) 



which can be viewed as a discrete analog of equations H4.3(l . This also shows that the difference 
vector Mfe_|_i — M^. is orthogonal to the rotation axis directed along (51, 0) G M"^, as expected. 
As a result, the map B : MP^ —> KP'^ given by H().13|l consists of the following 3 steps: 



1) . Given original set qi,q2,qo = ■\/ 1 — — i oi^^ finds components of Mfe from 16.1211 and 

of nlMkVLk from (|HT^ . 

2) . Given the components 



one finds new qi , (72 by solving the system of two algebraic equations originating from 



(Mfe+i)i = ((J22 + Jsa)?! - Ji2q2) y'^-qf-qj- (Jisqi + J23q2)q2, 



fe+l )2 



(6.16) 



((Jii + J33)'72 - -^1291) y'^-qf-q2 + (Jisqi + J23q2)qi- 



In R'^ = (91,(72,90) these equations describe two centrally symmetric quadratic surfaces 
Qi,Q2 which intersect the unit sphere qf + q2 + qi ~ ^ along curves Ci, C2 respectively. 
Each curve is a union of two ovals, which are centrally symmetric to each other. The 
intersection of Ci , C2 gives 4 complex points and 2 or none real points on P^ . Thus there 
are at most two different real solutions {qi\q2\qo^) with > 0. 

3). One chooses a solution (^q[^\ q2^\ q^^'' > 0^ and finally finds the last component {Mk+i)3 
by the formula 

(Mfe+i)3 = (Jii - J22)qiq2 - Ji2{ql - ql) - {Juqi + J23q2)qo, 

which is obtained from H6.16(l by substitutions k ^ k + 1 and the <z — > g. 



20 




Figure 6.2: A Circular Section of the Momentum Surface 



As a result, for n — 3 the map Mk —* M^+i given by (|6.13l) is generally ^-complex valued and 
2-real valued. In order to choose one of the 2 real branches, we must use some extra arguments, 
like existence of an additional integral, or, at least, to restrict ourselves with sufficiently small 
9i, 92, which correspond to rotations by a small angle 9. In this case only one of the solutions 
{qi\q2^) will be small and it is natural to choose it. 

It appears that the constrained energy integral H4.6|l of the continuous Suslov system is 
preserved by the discrete system as well. 

Theorem 6.2. The discrete Suslov system f6.1ci\) has quadratic integral 

{.hi + .h3)Ml + 2Ji2AfiA/2 + (J22 + J33)M|, (6.17) 
which gives rise to the following quartic integral in terms oj the parameters go, 91, 92.' 
H = (( J22 + Jzz)ql - 2Ji2qiq2 + (Jll + Jss)?!) 

• ((J1391 + J2392)' + [(Jll + J33)(J22 + J33) - J?2]9o) ■ (6.18) 

The proof is straightforward: substituting expressions (|6.12() and (|6.14() into H6.17|l gives the 
same expression in terms of go , 91 , 92 . 

The fact that (|6.17l) does not depend on M3 is quite natural: different branches of the map 
H6.13|l have the same value of the integral. 

It should be emphasized that the complete energy integral (Af, I^^Af) of the continuous 
Suslov problem is not preserved in the discrete setting. 

Invariant curves. As follows from Theorem 16.21 the map has invariant curves, which are 
either intersections of the sphere {qf + g^ + 9o = 1} with a quartic surface H{q) = h or, in the 
momentum space so* (3), intersections of the generalized quartic Steiner surface U with elliptic 
cylinders defined by H6.17|l . Thus, the invariant varieties are algebraic curves of order 8. 

Assume that quadratic form (J22 + ^33)91 ~ 2J129192 + (^33 + Jii)q2 is positive definite. 
Then, as follows from H6.18|l . on the upper hemisphere < 90 < 1 real invariant curves consist 
of two branches: for small positive values of h one branch is a small oval around the origin (0, 0) 
whereas the other branch is an oval close to the equator {go — 0} of the sphere. It may or may 
not intersect the equator. In the first case the opposite points of intersection are identified. 
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These different branches correspond to the two connected components of the intersection of 
the Steiner surface lA with the cyhnder. 

As value of the integral increases, the branches approach each other: the smaller one becomes 
bigger and the bigger shrinks. At a certain critical value h = h* the branches intersect at two 
opposite saddle points and form a separatrix, and for the next critical value h** > h* the two 
branches shrink to opposite center points. There are no real invariant curves for h > h** . Note 
that for h — h* and h = h** the elliptic cylinder is tangent to the surface U. An example of the 
invariant curves foliation is given in Figure 




Figure 6.3: Invariant Curves and the Equilibria Line on RP^. 



Remark. As noticed in ^Hl) in the absence of nonholonomic constraints, the map Mk Mk+i 
given by the discrete Euler-Poincare equations (|6.13|) is multi-valued, because, in general, the 
equation Mk = ^kJ — J^'k more than one solution. 

In presence of the constraints H6.10|l . the latter equation has generally a unique solution 
(except the points on self-intersection on U), however, as we saw above, the choice of A^+i or 
(Mj,+i)3 is not unique, and the map describing the discrete Suslov problem is multi- valued as 
well. 

Stationary solutions of the discrete Suslov problem. As follows from (|6.15() . if the 

initial values qi, q2 satisfy the condition Jisqi + J23q2 = 0, then 

(Mk+ih - (Mfc)i, (Affc+i)2 = (Mfc)2, {^iMk^kh = -{Mlh. 

that is, the coadjoint action i-^ VlJ^Mk^k is the mirror reflection with respect to the plane 
M3 = 0. Then it is natural to choose the multiplier Afe such that {Mk+i)z = {Mk+i)^- 

As a result, one of the branches of the map B has a one-parametric family of stationary 
solutions (equilibria) characterized by points of the line 

P = {^2 ^ { Jl3<Zl + J23g2 = 0}}/!?. 

They correspond to discrete versions of permanent rotations of the body in the classical Suslov 
problem. (In Figure Wi^ the set of equilibria points is represented by a straight line segment.) 
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In view of Ht).10|l . opposite points (91, 52, 9o) and (—(71, —(72, 9o) on P correspond to mutually 
inverse finite rotations and rf" respectively. 

As also follows from (|6.15l) . there are no equilibria points outside of this line. In particular, 
neither the saddle points nor the centers of the invariant foliation on MP^ are stationary points. 

Finally, note that, like in the continuous system, for a balanced inertia tensor J13 = J23 — 
all the solutions of (|6.15|) are stationary, i.e., the discrete body momentum Mk is preserved. 

Remark. The foliation of MP^ by invariant curves gives us a natural way of choosing the 
branches of the map B in the general case. Namely, if the initial point {qi, (72) lies in the domain 
S C KP^ defined by the condition < h < h* , i.e., it represents either a relatively small or 
sufficiently big finite rotation fi, then the points (51, (72) and (gi, ^2) have to belong to the same 
connected component of the invariant curve. In other words, if the initial point lies in the 
interior (exterior) part of S, one has to choose a real solution of Ht).16|l that has the smallest 
(largest) norm (7^ + ^21 respectively. 

On the other hand, if {qi, (72) lies in complement MP^ \ S, i.e., it is between the separatrices, 
then a real initial point (91,(72) may lead to complex ((71,^2) only. In particular, when the 
initial point is a center, the next point is necessarily complex, although the value of the integral 
remains to be real. 

If branches of the map EP^ 1-^ are chosen according to the above way, then the discrete 
time dynamics inherits all the main properties of the continuous Suslov problem. 

Namely, let A_ and A+ denote semi-planes of KP^ defined by conditions Juqi + J23'72 < 
(respectively > 0) and let 8_ and be semi-planes given by 

(J12 Jia + J22J23 + '^23'^3)'7i ~ {J11J13 + J12J23 + •^13^33)92 < 0, 
respectively > 0. 

Theorem 6.3. // the initial point q = (91,(72) Uss in the interior part of S <Z MP^, then for 
k —00 and k +00 the sequence {qfe} remains on the same branch of invariant curve 
and tends to the unstable equilibria semi-line P„ = P n 0_ and the stable equilibria semi-line 
Ps = P n 0- respectively. It lies entirely in one of the semi-planes A±. 

For the foliation indicated in Figure 16.31 the corresponding discrete time dynamics in the 
neighborhood of the origin is given in Figure where stable and unstable equilibria points on 
P as depicted as dots and circles respectively. 



0.2 - 



0.1 







-0.1 



-0.2 




-0.3 -0.2 -0.1 0.1 0.2 0.3 



Figure 6.4: Discrete Dynamics near tlie origin of RP^. 
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As follows from Theorem 16.31 for k —>■ — oo and k +00 the limit finite rotations Qk are 
mutually inverse. This property gives a perfect discrete analog of limit permanent rotations in 
the classical Suslov problem. 

Proof of Theorem \().l^\ First, we describe the discrete dynamics on the part Uq of the momentum 
surface U bounded by the condition 

E = (Jii + J33)Mf + 2J12M1M2 + (J22 + J33)M| < h*. 

For this purpose introduce a new coordinate system 

Ml = (Jl3(Jll + J33) + J12 J23)Mi + (J23(J22 + J33) + ^12^13)^2, 
M2 = J23M1 - J13M2. 

In view of relations (|6.12(l one has 

Ml = {Jmi + J23q2) [Q + {{Jii + J33){J22 + J33) - ^12)90] , (6.19) 
M2 - Qqo - iJi3qi + ^2392)', (6.20) 

Q = {J12J13 + J22J23 + J23J33)qi — {J11J13 + J12J23 + •^13^33)92- 

Using the properties ( Jn + J33)( J22 + >/33) ~ J12 > 0, 90 > 0, one can show that in the domain Uq 
the expression in square brackets in (|6.19l) is positive. Hence, on the segment of the line A^i = 
in 2Aq on has Ji3gi + J2392 = and it consists of stationary points of the map. The points of 
Uq with positive (negative) M 1 correspond to the points on the interior part of 5 C KP^ with 
positive (respectively negative) values of Jis^i + J2392. Next, in view of H6.15|l . 

M2,k+l = M2,k + (JuQl + J23q2f, 

which implies that the coordinate M2 always increases while the point approaches the line 
Ml = Q along the ellipse E{Mi, M2) = const. Then, as follows from (|6.20l) . for k —00, one 
has M2 < 0, Q < and for k —> 00, M2 > 0, Q > 0. As a consequence, the equilibria positions 
on P„ = P n 0_ are unstable and those on = P H 0+ are stable. 

Further, due to H6.15() . Mi,k+i — Mim = — («/i39i + ^2392)^ and, therefore, 

Ml,k+1 = (Jl3i?l + J23'?2)((Jll + J33){J22 + J33) - Ji2)<lQ- 

The latter and Ht).19|l implies that, unless Jisqi + 0/2392 = 0, the coordinates Mi.k and 
A4i^k+i always have the same sign, i.e., the sequence {Mk} lies entirely in one of the domains 
Uq n {A^i ^ 0}. Reformulating these properties for the interior part of the domain S C KP^, 
we arrive at the statement of the theorem. 



7 Discrete Unbalanced Chaplygin Sleigh 

Now we pass to discretization of the EPS equations (|5.3|) on the coalgebra se*(2). 
The two subsequent positions of the sleigh are given by the matrices 

(cos 9k -sin^fc Xk \ / cos^fe+i -sin 61^+1 Xk+i \ 

sin9k cos 9k yk , Xk+i = sin 6*^+1 cos^fc+i yk+i 
1/ \ 1 / 

The helical displacement in the body frame is defined by (2k — ^k^^k-\-i G SE{2) and straight- 
forward computation shows that 

(cos(A6'fc) -sin(A6'fc) cos 6*4; Axfe + sin 6*^ Ay^ \ 
sin(A6'fe) cos(A0fe) - sin 6'fc Ax^ + cos 6*4; Ay^ , (7.1) 
I ) 

A9k ^ 9k+i - cos 9k, Axk^Xk+i~Xk, Ayk =^ yk+i - Vk- 
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Following the expression 1)3. define the left-invariant discrete Lagrangian on SE{2) x SE{2) 
by replacing the helical velocity f in H5.8|l with X^^{Xk+i — Xk). Up to an additive constant, 
we get 



Ld(Xk+i,Xk) = itr {f2kSn^) - itr (jr2^ + 



(7.2) 



J/2 + ma^ mab 



mab 
ma 



J2 + mb"^ mb | , 
mb m 



where, as above, a, b are coordinates of the mass center C in the body frame and J is its scalar 
moment of inertia with respect to the origin A. This yields the following scalar expression 

Ld =^/^yl + yAxfe + (j + ma^ + mb^) (1 -cosA6'fe) 

+ am[{s\Ti6k+i - sm9k)Ayk + {cos9k+i - cos9k)Axk] 

+ bm[(cos0k+i -cos6k)Ayk - (sin 6*^+1 - sin 6'fe)Aa;fe]. (7.3) 

In the continuous limit, when 

A9k^euj + 0{e^), Axk ^ ex + 0{e^), Ayk = ey + 0{e^), e«l, (7.4) 
cos6'fc+i — cos^fe — — ecjsin6' + O(e^), sin^^^^+i — sin6'fc = etjcos^ + 0(£^), 

expression (|7.3|l divided by s transforms to the continuous Lagrangian l|5.1|l plus higher order 
terms in e. 

According to definition H3.2|l . the discrete momentum in the body 
Pk = {pe,k,Pi,k,P2,k) e se*(2), has the form 



Pl,k 



d_ 
Ik 
d_ 

d 



Ld{Ok + e,0k+i,Xk,Xk+i,yk,yk+i) 



6=0 



Ld{0k,Ok+i,Xk + ecos6k,Xk+i,yk + esin9k,yk+i) 



that is, 



P2.k = --^LdiOk,Ok+i,Xk - esm9k,xk+i,yk + e cos9k,yk+i) 



Pe,k = {J + ma^ + mb^) sin(A6'fc) + amV2,fc — bmVi^k , 
Pi,k = mVi^k - am{l - cos(A6'fc)) - 6to sin(A6'fc), 
P2,k = mV2,k + amsin{A9k) ~ bm{l - cos{A9k)), 



e=0 



e=0 



(7.5) 



where 



Vi,fc = (J?/c)i3 = AxkC0s9k + Ayk sin9k, 
V2,k = (J^/c)23 = -Axk sin 9k + Ayk cos 9k 

are "discrete velocities" in the body frame. 

Next, the coadjoint action on se*(2) can be written in form 

Piu.k ~ P2,kVi^k +Pl,fc^2,fe 

Adli^Pk = \ cos{A9k)pi,k + sin{A9k)p2.k 
~ sm{A9k)pi^k + cos{A9k)p2,k, 



(7.6) 



(7.7) 



In the absence of constraints the dynamics of the 2-dimensional body can be represented by 
the discrete Euler-Poincare equations 



Pk+i - Ad^Pfc, 



(7.8) 
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which gives the momentum conservation law written in the body frame. In particular, for 
a — b — (the mass center C lies at the origin), the system H7.8() . 1)7. 7|l yields 

sm{ek+i - 6k) = sin(6'fc - 6u-i), 
Axfc+i cos6'fc+i + Ai/fe+i sin 6*^+1 = Ax^ cos6'fc+i + Ay^ sin 6*^+1, 
-Axfe+i sin 6*^+1 + Aj/fc+i cos6'fe+i = -Axfc sin6'fc+i + Ay^ cos6'fc+i, 

which implies that for small 6 's the differences O^+i — Ok, Xk+i — Xk, and yk+i — Vk are the same 
for any integer fc, the result one expects from studying the continuous problem. 

Discrete constraint on SE{2). We now impose discrete left-invariant constraints on 
SE(2) X SE{2) in the form of restrictions on discrete helical velocities f?fc — Xk+iX'[ . By 
analogy with continuous constraint defined by (|5.11|l , a naive choice of a discrete constraint is 
just to set 

(r2fc)23 = - sin Ok Axk + cos 9k Ayk = 0. (7.9) 

This choice however is not the right one. Indeed, following our approach to discrete left-invariant 
constraints, admissible rotations and translations must be exponents of the matrices of the form 
H5.1()|l . In this case t) generates the subgroup SE{n — 1) and, according to ProDOsition l3.1l expO 
must be a covering of the homogeneous space SE{n) / SE{n — 1). 
In the particular case n = 2, when Qk is given by H7.1|) . we have 

Proposition 7.1. The variety S — expO C SE(2) is dijfeomorphic to the the canonical line 
bundle tt : £ — > — (zi : Z2) (Moebius cylinder) such that 'k~^{zi : Z2) — {v{zi, Z2),v G M} 
and it is defined by the condition 

[223 1 ~ [^11 



The latter yields the following constraint 



(7.10) 



— {Axk cosOk + Ayk sin^^) sin(A6'fc/2) 

+ (-Axfc sin 6ife -HA?/fc cos 6'fe)cos(A6'fc/2) = (7.11) 

or, equivalently, 

Fi,fc[l - cos(A0fc)] - V2M sin(A0fe) = 0. (7.12) 
The corresponding left-invariant constraint on SE{2) x SE{2) has the form 

- sm (xk+i - Xk) + cos (yk+i - yk) = 0. (7.13) 



V 2 

Observe that in the continuous limit (|7.4() this yields the constraint — isin6' -I- ycos9 — 
Proof of Provosition U . II For an element 5 G we have 

(0 —uj V \ / cos Ljt — sinwt ^sm{ujt) 

uj , exp(S't) = sinLot cosujt ^{1 — cos{ujt)) 
000/ \o o'^i 

where uj,v are arbitrary. As a result, for the points of the admissible shift subvariety, relation 
H7.10|l holds. Next, since 

, .a A ^23 1-cosujt AOk . . 

ujt = AOk, and — — ee ■ — = tan — -, (7.14) 

j/13 smwr 2 

in view of (jZU, we have lf7TI|l and (fTTHl . 

Finally, as seen from the last relation, the angle A^a; determines the quotient i723/r2i3, i.e., 
a line in — (i723, ■'^13) • As A^^ changes by 27r, the line rotates by tt, hence S is diffeomorphic 
to the Moebius cylinder. 
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Remark. As seen from relation H7.11|l . the matrices from S C SE{2) describe "circular 
translations" of the sleigh along the axis Xk of the blade: the points {xk, Uk) and {xk+i, Vk+i) in 
must lie on a circle such that the lines Xk and Xk+i are tangent to this circle. This property 
also implies that 

Axk cos 9k + Ayk sin dk = Axk cos 61^+1 + Ajjk sin 9k+i , (7 15) 

-Axk sin 9k + Ayk cos 9k = Axk sin 9k+i - Ayk cos 9k+i 

(see Figure inj. 




Figure 7.1: The geometry of the incremental displacements for the Chaplygin sleigh. 



The above constraint has also the following interpretation: in order to transfer the sleigh 
from {9k,Xk,yk) & SE{2) to {9k+i,Xk+iTyk+i) & SE{2) (assuming that this transition is pos- 
sible), one needs first to perform the rotation over A9k/2 at {xk,yk), which aims the sleigh 
towards {xk+i,yk+i), then slide the sleigh from {xk,yk) to (xfc+i, j/fc+i), and then perform an- 
other rotation over A9k/2 (now at {xk+iTyk+i))- 

Note that under the constraint (|7.12|) the image of the discrete Legendre transformation 
H7.5|l is an algebraic quartic subvariety U in se*(2) = (pe,pi,p2) and to a generic pair {pg,pi) 
there correspond four distinct points on U and four inverse images on 5 C SE{2). 

Discrete momentum locus 14 C se*{2). Below we concentrate on the important case 
b — 0, when the structure of the real surface U C se*(2) becomes simpler. It is more convenient 
to describe the image U oiU inM.^ ^ {pg,pi, z), where pi = api + 2ma^, z — sin(A6'). 

Lemma 7.2. 1). For b — the surface hi is given by cubic equation 

n{pe,Pi,z) = J2z3 _ 2Jpe,kZ^ + [pi + 2Jpi+pl)z - 2pgpi = 0. (7.16) 

U lies entirely between the planes z — ±1 and is tangent to them along the lines £± = 
{±P0 — pi = J} respectively. The pg- and pi- axis belong entirely to hi. 

2). For the parts of hi over the quadrants 

L++ = {-pi + pg > J} n {-pi - pg > J} and 
L— = {-pi +pg < J}n{-pi - pg < J} 

one has cos(A0) > 0, i.e., —tt/2 < A9 < 7r/2 and in the rest of quadrants one has 
cos(A6l) < (7r/2 < < 37r/2). 
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3). The projection H : [/ ^ = {ve,Pi) is one-to-one except the above segments and the 
interior of triangular domain bounded by the discriminant curve 

p\ + &JpI+pI{12J^ + 2pl)-p^{lOJpl - 8J3) _ j^pI = q. 

The curve is symmetric with respect to pi-axis, it is tangent to pg-axis at the origin (0, 0) 
and has 3 cusp points with coordinates (0,-2 J), (ci, C2), (— ci, C2), with some positive 
constants Ci,C2. In this domain the projection W is 3 to 1. 
4-). The curve {Vi = 0} dU is projected onto the ellipse 

£ = {pe = (J + ma'^)sm{A0), pi = ma^ {1 + cos{A9)) \ A0 e {0; 27r)} . (7.17) 

Inside the ellipse the values of Vi are negative and outside are positive. 

Note that the pomt O with coordinates pg — 0, pi — 2ma? corresponds to the origin in the 
iPe,Pi) phase plane and in a neighborhood of this point the projection 11 is one-to-one. An 
example of the surface U for J = 1.5 is presented in Figure FT^ 




Figure 7.2: The surface U. 

Proof of Lemma [7.21 1). Using the condition H7.12|l . we exclude Vi,V2 from the first two 
equations of (|7.5|) to obtain the following condition on Adk, 

J sin^iAek) - pg,k sm{A0k) + {2ma'^ + api^k)[l - cos{A9k)] =0. (7.18) 

This equation always has trivial solution AO^ = 2mT, n e Z. Setting z — sin(A6'fe), cos(A6'fe) = 
z^, Pi = apx + 2ma^, we arrive at a quartic polynomial equation with respect to z, which 
has the root z = 0. Factoring it out and omitting the index k, one gets the cubic equation 

msi). 

Now setting in H7.16|l z = ±1, we get {J ^ pg -l-pi)^ = 0, which implies that U is indeed 
tangent to the planes z = ±1 along the lines £±. Finally, setting z — pg — or z = pi = 0, one 
sees that equation H7.15|l is satisfied for any pi and pg respectively. 

2). For fixed pg,pi, each root of H7.16|l gives a solution of (|7.18|) with a sign of cos(A6'fc) 
appropriately chosen. As seen from (|7.18|) . for large \px\ and small z = sin(A6'), the value of 
cos(A0) must be close to 1, whereas for large pg > and small \px\, cos(A0) must be negative. 
Since the sign of cos(A0) can change only under passage from one quadrant on the plane {pg, Y) 
to another one, this proves item 2). 

Items 3), 4) are verified by direct calculations. 
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Discrete dynamics on se*(2) with the constraint. According to (|3.8|) . the discrete 
Euler-Poincare-Suslov equations associated with the constraint (|5.11l) have the form 

Pfc+i = Ad^,^Ffc + Afc(0,0,lf . (7.19) 

Substituting here expressions H7.7|l . we find that under the constraint H7.12|l the first two 
components of Pk+i have the form 

Pe,k+i = (J + + rnb^) sin(A6'fc) — bmVi,k 

+ am[- Ax k sin 9 k+i + Aj/fe cos 6'fc+i], 
Pi,fc+i = mVx.k - 6msin(A6'fe) 

+ am[Axk cos ^fe+i + Ayk sin 6*^+1] , 
which, in view of H7.15|l . (|7.5() . yields 

Pe,k+i = POM - '2amV2.k , ^g) 
Pi,fc+i = Pi.k + 2am{l - cos{A9k)) . 
Expressions l|7.20|l . H7.12|l define a multi- valued map U oi S ^ S which consists of 3 steps: 

1) . Given AOk, Vi^k, one finds V2,fc form the constraint 1)7. 12f) and then pe,k,Pi,k,P2,k from the 

Legendre transformation (|7.5() . 

2) . One finds pe,fc+i,P2,fc+i from 1)7.20)1 . 

3) . One finds Adk+i, Vi^k+i by choosing a solution of the system of equations 

PeM+i = {J + ma^ + nilp-) sin(A0fe+i) + (am- — .^°^!^^^^X^^ - bni] Vx,k+i , 



sin(A6'fc+i) 

Pi,fc+i = mVx,k+i - am(l - cos(A6'fc+i)) - 6to sin(A6'fc+i), 
which are obtained from 17.5)1 . 1)7.12)1 by replacing k k + 1. 

Theorem 7.3. Equations fTW^ preserve the quantity 

E^mpl + 2bmpgpi + ( J + m{a^ + b^))pl, (7.21) 

which coincides with the truncated energy integral 1^5.^) of the continuous Chaplygin sleigh. 

Proof. Substituting expressions ()7.5)) and 1)7.20)1 into 1)7.21)1 and taking into account the con- 
straint 1)7.12)1 . one obtains the same expression in terms of AOk, Vi.k and V2.k- 

Since the quadratic form 1)7.21)1 is positive definite, we conclude that the invariant manifolds 
of the map 1)7.20)1 are the ellipses in the pepi-plane. 

Stationary solutions of the discrete Chaplygin sleigh. As follows from ()7.2U)) . for 

the initial conditions {Adk ~ 0, V2.A; = 0} one has 

Pe,k+i^Pe,k, Pi,k+i ^ Pi,k- 

Hence, it is natural to choose such Afe in ()7.19)) that P2,fe+i ~ P2,k as well. Thus, like the 
continuous system ()5.4)l . for a 7^ the map 1)7.19)1 has a family of stationary solutions which, on 
the momenta plane (pe,pi), is represented by the line {pe + hpi — 0}. Such solutions correspond 
to shifts in the [x, y)-plane along the axis of the blade by constant distances. 

On the other hand, for a = all the solutions are stationary. That is, in contrast to the 
case of absence of constraints, when the discrete momentum in space is preserved, now it is the 
momentum in the body P, which is preserved. In view of 17.5)) . this implies 

A6k+i — A6k, Vi^k+i — Vi,k- 

As a result, the discrete trajectory on the plane (a;, y) consists of displacements along a circle 
with radius p = Vi,k/ sin(A6'fe) ^ . The same behavior occurs to the continuous sleigh for a — 0. 



^As numerical simulations show, if one chooses the naive constraint 17.911 instead of 17.101 1. then for a — the 
discrete trajectory on the plane {x,y) lies on a spiral. 
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The case b = 0, a ^ 0. In this case the map {p0,k,Px,k) — > {pe,k+i,Pi.k+i) has a Hne 
of stationary points pg — 0, and, according to Theorem 17.31 the discrete trajectories he on 
symmetric invariant eUipses mpg + ( J + •ma?)p\ — E. Without loss of generahty, we assume 
a > 0. Then the following property holds. 

Theorem 7.4. In the neighborhood of the origin O hounded by the condition E < rn^ a? {J +ma'^) 
the map is single-valued and has a hi- asymptotic behavior similar to that of the continuous 
Chaplygin sleigh system. Namely, for k —oo, the point {pe^k,Px,k) approaches, along the 
corresponding invariant ellipse, a point of the segment {pg = 0, —ma < pi < 0} of unsta- 
ble stationary points, and for k — > +00, the point (pe.k^Px.k) approaches one of the points of 
the segment {pg = 0, < pi < ma} of stable stationary points. In both cases the sequence 
{{pe,k,Px.k)} remains in one of the half-planes pg < or pg > 0. 

Proof. Part (3) of Lemma lT. 21 implies that the map is single- valued in the region E < m'^a'^{J + 
ma^). 

Next, as follows from the first relation in H7.20|l for a > 0, the increment pi^k+i ~ Pi.k is 
always greater than or equal to zero. Then, to prove the bi-asymptotic behavior, it remains 
only to show that the sequence {{pg,k,Pi,k)} lies entirely in one of the half-planes pg ^ 0. 

Indeed, let the point {pg.k,Px.k) be inside the ellipse £ given by (|7.17|l . First, assume that 
Pe,k > 0. Then, in view of items (2), (4) of Lemma lT?^ and the constraint H7.11|l . Vi^k and V2.k 
are negative. According to H7.20|l . the increment pg.k+i ^ Pe,k is then positive. Similarly, for 
pg^k < one has pg,k+i - pg^ < 0. 

Next, if < pi,k < rna and {pg,k,Px.k) lies in the domain E < {J + ma'^)m'^a^ , then, using 
1)7. 5(1 . H7.11|l . one shows that 2amV2,k > Pe,k for pg^k > and 2amV2,fc < pg^k for pg^k < 0. 
Therefore, in view of (|5.3|) . pg,k+i > 0, respectively, < 0. 

As a result, in any case, pg^k and pg.k+i cannot have different signs, which proves the theorem. 

One concludes that in the neighborhood of the origin O the discrete-time dynamics is similar 
to that of the Suslov problem illustrated in Figure 16.41 

We conclude this section with an example of the discrete sleigh trajectory on the plane (x, y) 
compared to a continuous trajectory for 6 = with a cusp, as presented in Figure [^31 




Figure 7.3: A typical discrete sleigh trajectory. 



8 Conclusions 

The discrete nonholonomic Suslov problem and the Chaplygin sleigh that we introduced in 
this paper properties of their corresponding continuous-time dynamical systems; in particular, 
they preserve the reduced constrained energy and, in the balanced case, the momentum. This 
enables one to obtain explicit solutions for the momentum dynamics of both discrete systems 
in terms of theta-functions and exponents. It is not currently clear if the complete solvability is 
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due to the (low) dimension of the systems and if it is possible to construct completely solvable 
discretizations of the multidimensional Suslov and Chaplygin problems. These issues will be 
addressed in a future publication. 

On the other hand, by modifying our approach, one can also consider discretizations of 
nonholonomic LR systems on Lie groups. For such systems, the Lagrangian is left-invariant 
while the constraint distribution is ri^/ii-invariant. The discrete dynamics of such systems, as 
well as the existence of their invariant measure, is currently being developed and will be exposed 
in a future publication. 
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